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We propose a theory based on simple physical arguments that describes a non equilibrium steady- 
state by a temperature-like parameter (an "effective temperature"). We show how one can predict 
the effective temperature as a function of the temperature of the environment for a specific case 
of non-equilibrium behavior: radiation amorphization. The main idea that we present is that 
the amorphization process is inherently connected to the dynamical arrest that a liquid undergoes 
when it transforms into a glass. We suggest that similar arguments may hold also for the effective 
temperature under plastic deformation. 



The idea that the state of a system in a non- 
equilibrium steady-state obeys a Gibbs-like distribution 
with a temperature-like parameter, called the "effective 
temperature" was studied extensively starting from the 
pioneering work of Edwards [1 . Edwards assumed that 
a granular material under external perturbation explores 
a set of equal probability states and therefore the micro- 
states are described by a stationary probability measure 
analogous to the Gibbs measure, but where a configura- 
tional temperature (that he called "compactivity" ) re- 
places the thermal bath temperature. Cugliandolo et 
al.[2] showed, for mean field models, that such an ef- 
fective temperature has other properties desired for a 
thermodynamic description of non-equilibrium states, for 
example it controls the direction of heat flow. Further- 
more, it was shown that fluctuation-dissipation theorems 
are violated in a very specific way: at short time scales 
the fluctuations are related to the susceptibilities in the 
usual way, described by the fluctuation dissipation the- 
orem involving the thermal bath temperature. At large 
time-scales the usual fluctuation dissipation theorem is 
not obeyed but satisfies an analogous expression that in- 
volves the effective temperature. In later work similar 
behavior was observed for sheared glasses at tempera- 
tures below the glass transition (about 66% T g ) [3] and 
in sheared foam at zero temperature This suggests 
that a statistical mechanics description with a Gibbs- 
like measure might apply for these systems. It was also 
found that the response of amorphous materials to high 
shear-rates is related to their structure by an exponen- 
tial function, which suggests that the structure might 
be described by a Gibbsian-like measure, similar to the 
ideas of Edwards [6l [7] . Two important theories of plas- 
tic deformation in amorphous solids, the Shear Trans- 
formation Zone theory (STZ) |5] and the Soft Glassy 
Rheology theory (SGR) |S} use the idea of an effective 
temperature as a control parameter. In STZ theory it is 
assumed that the distribution of soft areas in a material 
obeys a Gibbs-like measure. Bouchbindcr et al.JTU] have 
extracted an effective temperature from non equilibrium 
values of an order parameter by inverting the equilibrium 



equation of state of the order parameter (assuming that 
the equation of state is still valid). In Boue et al. [TT] 
this approach was taken one step further by showing ex- 
plicitly that several coarse-grained structural properties 
of at least two different models of amorphous solids un- 
der extreme shear (plastic deformation) are described by 
a Gibbsian measure controlled by an effective tempera- 
ture. The value of the effective temperature as a function 
of the temperature was measured and was found to show 
non-monotonous behavior. It was also shown that the 
stress and the potential energy depend on the effective 
temperature. It is usually not possible to predict the ef- 
fective temperature theoretically, which leads to the use 
of fitting functions or scaling arguments. 

In this work we show, using molecular dynamics simu- 
lations, that the steady-state structure of an amorphous 
solid under irradiation can be described by a Gibbs-like 
distribution. More importantly, we show how one can 
understand and predict the value of the effective temper- 
ature at the steady-state by connecting the fast and slow 
dynamics in the system in a way that is deeply rooted 
in the nature of the glass transition. Our explanations 
and theory are different than the ones given by Martin 
in [12] , who predicted an infinite effective temperature at 
zero thermal bath temperature and also from the argu- 
ments given by Lund et al. in [T3] who suggested that the 
effective temperature has a finite value at zero bath tem- 
perature and from Bouchbinder et al. [H] who suggested 
a thermodynamic explanation for the non-monotonic be- 
havior of the effective temperature as a function of the 
thermal bath temperature. 

The structure and slow dynamics of a super- 
cooled liquid: In a series of recent papers (cf. [T51H6] ) it 
was shown that the structural changes that a supercooled 
liquid undergoes on cooling can be represented by a fi- 
nite set of species which are indexed by 1, 2, • • • , n. The 
precise nature of these species may change from model to 
model, but they are always formed by particles and their 
nearest neighbors. The main advantage of these species 
is that they obey discrete statistical mechanics, in the 
sense that their temperature-dependent concentrations 
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(Cj)(T), are determined by a set of constant degenera- 
cies gi and enthalpies Hi such that: 
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After applying this description we always find that some 
species vanish at temperatures that are close to the glass 
transition temperature. This qualitative observation was 
made quantitative by summing the values of one or more 
of the concentrations of the vanishing species to what is 
called the 'liquid- like' concentration (Ce)(T). The inverse 
of this concentration provides a length-scale (the typical 
distance between 'fluid' species): 

f(T) = [(C t ){T)]- 1/d , £ ^ °o when T ^ 0, (2) 

where d is the space dimension. It was amply demon- 
strated on a large variety of models that the relaxation 
time T a (T) measured using correlation functions in the 
super-cooled regime is determined by this diverging scale 
according to 
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where fx is a typical free energy per particle and r 
a typical vibration time. In the following we will use 
the results of this analysis for a specific model: the 
ST (Shintani-Tanaka) two-dimensional liquid model |17j . 
The ST model has N identical particles of mass m; each 
of the particles carries a unit vector that can rotate on 
the unit circle. The particles interact via the potential 
Ufa,9 h 9j) = Ufa)+AUfa,6i, 3 ). Here Ufa) is the 
standard isotropic Lennard- Jones 12-6 potential, whereas 
the anisotropic part AU(nj,9i,6j) is chosen such as to 
favor local organization of the unit vectors of two in- 
teracting particles, such that they form 126° with the 
radius-vector between them. This kind of interaction 
prefers structures with a five-fold symmetry, which frus- 
trate crystallization. For full details of this model the 
reader is referred to Refs. [I7Hl9| : here it suffices to 
know that with the parameters chosen in Ref. [T7j the 
model crystallizes upon cooling for A < 0.6 whereas for 
larger values of A the model exhibits all the standard 
features of the glass transition, including a spectacular 
slowing down of the decay of the correlation functions of 
the unit vectors C R (t) = (1/iV) Ei( u * W ' u i(0)), which 
is very well described by Eq. In [HI US] and pi] the 
above statistical mechanics theory was constructed for 
the ST model. We chose to separate the particles into 
three species, according to the possible number of neigh- 
bors interacting in the most favorable way with a given 
particle: particles with two favorable interactions are des- 
ignated as belonging to the "blue" species (n = 2), par- 
ticles with one favorable interaction are labeled "green" 
(n = 1) and ones with no favorable interactions at all are 
labeled "red" (n = 0). In previous work we have found 



the exact numerical values of the constants gi and Hi and 
showed that equation [I] describes accurately the concen- 
trations of species measured from experiments, and that 
equation [3] describes the typical relaxation times of the 
autocorrelation function Cn{t) [T71 [T5] - 
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FIG. 1: (Color online): Representative configuration of par- 
ticles in the ST model with examples for the three different 
species. Particle with two favorable interactions are assigned 
as a "blue" species, particles with one favorable interaction 
are labeled "green" and ones with no favorable interactions 
at all are labeled "red" . 



Simulations: The molecular dynamics protocol we 
implemented was the following: First, we carefully equi- 
librated a large number of independent configurations 
with N — 1024 particles in the NVT ensemble using the 
Berendsen thermostat over a wide range of temperatures. 
After this, we turned our equilibrium super-cooled liq- 
uids into amorphous solids by minimizing their poten- 
tial energy using the conjugate gradient algorithm. This 
procedure can be thought of as quenching a liquid in- 
finitely fast into a disordered solid whose temperature 
is formally T = 0. Then we irradiate at a constant rate 
while maintaining constant thermal bath temperature Tf,. 
We "irradiate" the material by giving a particle a very 
high momentum at random. This simulates the effect 
of an energetic particle bombarding the material. The 
high momentum creates a cascade event where the parti- 
cle transfers its energy to the surrounding particles. The 
simulations were performed at a constant temperature 
with a constant rate (average time between hits) and in- 
tensity (average momentum transfer) of radiation. The 
rate of bombardment of the material was much slower 
than the rate of thermal fluctuations and was chosen such 
that the material would not heat and melt. 

Effective temperature: In [TT] we found that the 
concentrations of species in steady-state plastic defor- 
mation (Ci) ss are described by the following Gibbsian 
distribution: 
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where gi and Hi have the same values as in the equi- 
librium Gibbs distribution and the effective temperature 
T e jf is a parameter that depends on the external forc- 
ing and thermal bath temperature in a non-trivial way. 
In figure [2] we show the concentrations (Ci) measured 
from simulations in an initial amorphous configuration 
(squares) and the concentrations measured at the steady- 
state under irradiation (circles). The curves are the the- 
oretical concentrations given by the equilibrium Gibbs 
measure for this system - equation [T] The concentra- 
tions in the initial quenched state follow the Boltzmann 
distribution. This is an indication that the structural 
aspects that are reflected in the species are not sensi- 
tive to the quench. Under irradiation the concentra- 
tions undergo a significant change (see inset) and reach 
new steady-state values. Remarkably, the new concen- 
trations, despite having changed significantly, still cor- 
respond to values of equations 1 with some "effective" 
temperature. For different radiation rates and thermal 
bath temperature the steady-state concentrations reach 
different values but they all agree well with the equilib- 
rium Gibbs distribution with some effective temperature. 
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FIG. 2: (Color online): Concentrations of species with two 
strong bonds' (blue), one strong bond (green) and no bonds 
at all(red) as measured from simulations against the equilib- 
rium theory (smooth lines, colors the same as concentrations). 
Squares represent the concentration at an initial random con- 
figurations quenched from a temperature T = 0.29. Circles 
represent the concentrations of the species at non-equilibrium 
steady-state under radiation. Inset: concentrations as func- 
tion of time starting from an amorphous solid. 



Theoretical description of the steady-state ef- 
fective temperature: In figure [3] we show the effective 
temperature measured at the steady-state for a given ra- 
diation rate as a function of the thermal bath tempera- 
ture (blue circles). In the inset we can see the steady- 
state values of the effective temperature of the same 
model under plastic deformation at two different, con- 
stant strain rates (see [llj). One can observe that in 



both cases the behavior is quite similar: At high tem- 
peratures the dynamics is dominated by the fluctuations 
imposed by the thermal bath and the irradiation has little 
effect. In that regime the equilibrium Gibbs distribution 
accurately describes the statistics. As the temperature is 
lowered, the rate of irradiation (strain-rate in the case of 
plasticity) becomes comparable to the time scale associ- 
ated with thermal fluctuations (eq. [3]). At that point the 
effective temperature starts to deviate from the thermal 
bath temperature. What happens at lower temperatures 
is completely un-intuitive: as the temperature is lowered, 
the effective temperature decreases moderately and than 
increases again until it reaches a finite effective temper- 
ature at zero temperature. But what is determining the 
values of the effective temperature? In the following we 
will give a physical explanation of this behavior. The 
basic idea underlying our analysis is that whenever an 
energetic particle transfers its energy to one of the par- 
ticles in the material, it transfers some of its momentum 
to the neighboring particles, which in-turn transfer some 
of their momentum to their neighbors and so on, in a 
"radiation cascade" . In many different materials, this is 
equivalent to a local melting of the material around the 
point of impact. Since the rest of the material is kept at 
a much lower temperature, the kinetic energy dissipates 
rapidly and the region cools down. Since the cooling pro- 
cess is very fast it will result in a local glass transition. 
We can describe this process using two rate equations 
(based on Fourier's cooling law): 



T K = - — [T K -T b ], 
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The first equation describes the relaxation of the kinetic 
temperature T K of the molten region. After long enough 
time, all of the kinetic energy will dissipate and the ki- 
netic temperature will equilibrate with the thermal bath 
temperature T b . The rate of relaxation is controlled by 
the thermal conductivity of the material and is mate- 
rial dependent. The second equation describes an aging 
process. At times longer than the relaxation time T a (T) 
the structure of the molten area - described by the ef- 
fective temperature T e tf - reaches thermal equilibrium 
(note that one can use instead of equation [3j the Vogel 
Fultcher expression to a similar effect). The solution to 
the first equation is: 



T K (t, T b ) = [T K (t = 0) - T b ]e-*/ T *" + T b 
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and to the second equation: 
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By choosing a cutoff time r c / we get an expression for 
the effective temperature as a function of the thermal 
bath temperature T eff (Tb) = T eff (r c f,Tb). This equa- 
tion is analytically solvable if we approximate and as- 
sume that T K = Tfc. However, for small temperatures 
this assumption fails and the approximation is not accu- 
rate. The non-Arrhcnius increase of the relaxation time 
at low temperatures means that the molten area will not 
reach the equilibrium structure but will undergo a glass 
transition and get "stuck" in some metastable state. In 
order to find the effective temperature that corresponds 
to this state, we numerically integrate the above equa- 
tions with a high initial kinetic temperature which we 
guess to be: T K (t = 0) = T eff (t = 0) = 0.5, repre- 
senting the temperature of the molten region just after 
a cascade event. We also chose r t h = 200 and introduce 
a time-cutoff t c / = 400 after which we assume that the 
structure does not evolve. The time scale r c / is related 
to the frequency of cascade events since this frequency 
controls the typical time between re-melting of a specific 
region. In figure [3] we observe the agreement between 
a numerical solution of the equations and the data from 
simulations. The behavior of the effective temperature as 
a function of the thermal bath temperature is captured 
by the theory. This seems to apply, at least qualitatively, 
to plastic deformation. Despite being a very different 
process, plasticity in amorphous solids shares similarities 
with radiation damage: plastic events tend to be local- 
ized, a large amount of kinetic energy is released in the 
process from these localized regions, and the events result 
in irreversible structural changes. 




FIG. 3: (Color online): Effective temperature measured from 
simulations with the same irradiation power but different 
thermal bath temperatures. The yellow line is the theoretical 
result. Inset: same result for simulations of plastic deforma- 
tion with the same model at two different strain-rates. 

Discussion: We have constructed a theory for the val- 



ues of the steady-state effective temperature of a material 
under irradiation. This result is based on a deep con- 
nection between the effective temperature and the glass 
transition which was proposed in the past |20j and is ex- 
plained here. Due to the obvious similarity to the behav- 
ior of the effective temperature in plastic deformation, 
we suggest that a modified theory will be able to predict 
the effective temperature in plastic deformation as well. 
Future work will focus on applying the same formalism 
to more realistic simulations and experiments and also 
on relating structural properties to the effective temper- 
ature, for example in irradiation induced flow (see for 
example [21]). 
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